Oscillatory vortex interaction in a gapless fermion superfluid with spin population 
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Using effective field theory approach we study a homogeneous superfluid state with a single 
(gapless) Fermi surface, recently suggested as a possible phase for an ultracold Fermi gas with 
spin-population imbalance. We find an unconventional form of the interaction between vortices. 
The presence of gapless fermions gives rise to an additional, predominantly attractive, potential 
oscillating in space, analogous to the RKKY magnetic interaction in metals. Our study then leads 
to an interesting question as to the nature of the vortex lattice in the presence of the competition 
between the usual repulsive logarithmic Coulomb and the fermion-induced attractive oscillatory 
interactions. 



Recently a beautiful array of vortices has been ob- 
served in the MIT experiment of ultracold Fermi gases, 
providing a definitive evidence for superfluidity in this 
class of systems. [l[ When fermionic excitations are ei- 
ther fully gapped or located at a few nodal points such 
as in a d-wave superconductor, the physics of vortices be- 
longs to the universality class of the XY model where the 
phase of the superfluid order parameter plays the domi- 
nant role. In this case, the vortex sector of the XY model 
is described by the logarithmic Coulomb interaction be- 
tween the vortices. 

The development of imbalanced (spin-polarized) cold 
Fermi gases is set to alter this standard picture in a fun- 
damental fashion, with a number of exotic superfluid 
states recently proposed or revisited 0, [3, I3- One of 
those states commonly found in various theoretical ap- 
proaches d, 0, 0] is a homogeneous gapless fermion su- 
perfluid with a single Fermi surface (FS) on the molec- 
ular side of the Feshbach resonance (the BEC regime). 
We will call this phase BPl (breached-pairing state with 
a single Fermi surface) after Ref. Q (also dubbed as 
"magnetized superfluid" and denoted SFm in Q). Two 
recent experiments on the polarized Fermi gases have 
brought the subject to the frontier in the cold atom 
physics . Several theoretical works have found an ana- 
logue of the BPl in the trap - the superfluid- normal- 
mixture phase. [13] It is then of great interest to examine 
properties of this novel superfluid. 

In the present work, we study the effect of gapless 
fermionic excitations on the vortices in the BPl phase 
by explicitly incorporating these excitations in the low- 
energy effective field theory. The resulting interaction 
between vortices is no longer of the pure Coulomb form, 
but contains an additional fermion-induced contribution 
that is attractive (between the vortices with the same 
vorticity) and oscillates on a length scale set by the spin 
polarization, analogous to the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) magnetic interaction in metals. 



The Bogoliubov-quasiparticle energy spectrum in a 
system of two fermion gases with an equal mass and un- 
equal chemical potentials is given by [2| {h = 1) 
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where fi = (/i| + /ix)/2 is the average chemical poten- 
tial and (5 = /i| — /i| the chemical potential mismatch. 
Since our treatment concerns the superfluid phase real- 
ized on the BEC side of the Feshbach resonance, /i < 
in what follows. For definiteness, we assume that S > 0. 
If S is sufficiently large (i.e., if (5/2 > y/Jj?~\~A^), the 
lower branch of the above dispersion is gapless, with 
a single effective FS. For brevity, we hereafter denote it 
as £k- When linearized in the vicinity of the effective 
FS, this gapless dispersion adopts the generic form Ek — 
Wb(|k| — fcfc), parameterized by the "Fermi velocity" Vh 
and the radius kh = [GTr^nf,]^/^ of the "breached-pairing 
Fermi ball" in momentum space, where rif, = — 
(here ni > n^, as follows from (5 > 0). 

In the present context, the superfluid phase field rep- 
resents the phase of the complex Cooper pair ampli- 
tude ((V'T'0i) = l('/'TV'i)|e*^) and can be separated into 
a (singular) vortex part and a smoothly fluctuating field 
(j) via V6'(x, r) = —a + V0(x, r), where a is the vor- 
tex gauge field. Let us consider a static configuration 
of vortex lines of unit winding number, all parallel to 
the z-axis, with a "vortex charge" density distribution 
p(x) = 27r^^(S(2'(x — Xq). The gauge field a is re- 
lated to p by the standard equation V x a = ~Kop{x.)ez, 
where kq ~ ■nh/m (h restored for clarity) is the circula- 
tion quantum. In momentum space this reads q x aq = 
iKop(q)e2, where q is a two-dimensional wave- vector 
(q • = 0) of the phase fluctuations and p{q) the 
Fourier transform of /o(x). For convenience, we adopt 
the Coulomb gauge V • a = in which the vector field aq 
is purely transverse (q • aq = 0). 

We start from a low-energy effective Lagrangian for 
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the gapless branch of the Bogohubov quasiparticles (de- 
scribed by the field xi^)) and the superfluid phase field. 
This Lagrangian obeys two global U{1) symmetries, one 
of which corresponds to the total atom number conserva- 
tion (J7c(l)), and the other one to the conservation of spin 
population imbalance (J7s(l)). Son and Stephanov Q 
first constructed an effective Lagrangian of this kind for 
the case of small gapless fermion density (corresponding 
to small spin population imbalance) to obtain a phase 
diagram. Here, we extend their approach to an arbi- 
trary spin population imbalance. (ll|In the imaginary- 
time path-integral formalism, our effective Lagrangian 
reads 



C= x*[dr + e{-iV)]x + ci{dref 
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where '• • • ' stands for higher-order derivative terms of 
the field; j = {x*^X - Vx*x)/(2mpi), where nip = 2m 
is the mass of the Cooper pair; e(— iV) is the gapless 
fermion dispersion, written in the coordinate representa- 
tion. The phenomenological parameters ci, C2, and C3 
are not constrained by the C/(l) symmetries, but can be 
determined from the Galilean-invariance of the fermion- 
independent part of the above Lagrangian ; ci = dn/dfi 
(where n = nij+ni), while C2 and C3 are related to the su- 
perfluid density in a manner to be described shortly. Un- 
der Galilean boost with velocity u we have the following 
transformation properties : dr d'^, = «(u-V), V 
V = V , x(x,r) ^ x'(x',r') = x(x, r) , 0(x,t) ^ 
6''(x',t') = 6'(x, t) — TTipU • X — ^nipU^T. An alternative 
effective theory, formulated in terms of majority-species 
original fermions rather than Bogoliubov quasiparticles, 
is given in Ref. 



Because the Lagrangian has to be invariant under the 
C/c(l) particle number symmetry 6 6 + it contains 
the coordinate and time derivatives of 9, but not 9 itself. 
Recast in terms of </> and a, it adopts the form 



C= X* [Or + e{-iV)]x + ci(9,(/.)2 + C2(V0 - a)^ 



+C3X*X 



idr4> ' 



(V0 - a)2 



2mB 



(V0-a).j 



(3) 



We derive an effective phase-only action S[9] = 
S[(f),SL] by integrating out the fermionic degrees of free- 
dom: e-'5M = J D{x*:X)e^'^^^''^^: where S[x..9] = 
dr J dx C is the Euclidean action corresponding to 
the last Lagrangian {(3 = {kBT)~-^). In order to ac- 
complish this, we first note that the fermion field enters 
the Lagrangian ^ through a quadratic form x*Kx = 
X*i-Qo' + X)x, where Go = [-Or - e{-iV)]-^ is the 
noninteracting fermion propagator, and X — X^^^ +X'-^^ 

where X^^) = cgi^r^ + t^(V0 - a)- V and X^^) = 
2^(V(/) — a)^ are respectively of the first and second or- 
der in fields and a. Integrating out the x fields gives 



rise to the fermion contribution to the action S'[0, a]: 
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SF\'i>A = -trlnif = const. + ^ -iY\((5^XY\ . (4) 
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We expand ^ to first order in X^^^ (tree level) and 
to second order in X'^^^ (one-loop order). The effective 
phase-only action 5'[(7!', a] is then obtained by gathering 
Sp\(^, a] and the fermion-independent terms of the origi- 
nal action. In momentum-frequency space (g = (q, zcj;)) 
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where the first and second terms correspond to the prop- 
agating Goldstone mode and the topological vortex part 
of the broken ?7c(l) symmetry, respectively. [Summation 
over repeated indices in the last equation is implicit.] 
Here Ii{q) is the fermion density polarization bubble, 
while 

^'^('?) = 47E^oW^o(^ + 9)(^» + f)(% + |X6) 
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represent the longitudinal and transverse current-current 
correlation functions, respectively, with k_L being the 
transverse component of the three-dimensional vector k 
with respect to q. 

It is important to point out that there is no RPA-type 
correction from the interaction vertex j • V(/> to the trans- 
verse current-current correlation function; this is mani- 
fest in our choice of the Coulomb gauge for the topolog- 
ical gauge field a. 

The magnitude of wave- vector q of the phase fiuctu- 
ations has an upper cutoff of order /ca — (2mA) a 
momentum scale corresponding to the pairing gap. This 
choice of momentum cutoff is a consequence of the BPl 
phase being realized in the strong-coupling regime on the 
BEC side of the Feshbach resonance, where A is related 
to the binding energy of a Feshbach molecule. It is also 
important to emphasize that typical |q| is not necessarily 
small with respect to kb\ the latter is controlled by the 
spin imbalance and can be arbitrarily small. 

It is straightforward to show that Rij{q) = R{q)Sij. 
Consequently, the phase-only action in Eq. (O in the 
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zero-temperature static limit reduces to 
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where Pq and i?q are the zero-temperature static hmits 
of P{q) and respectively. 

The superfluid mass density Ps ~ rigmp, which plays 
the role of rigidity in the present problem ("spin- wave" 
stiffness in the XY-model terminology), can be identi- 
fied from the long-wavelength (q — > 0) limit through the 
relation 
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This equation constraints the coefficients C2 and C3, which 
reduces the phase-only action in Eq. ([5]) to 
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which is free of the phcnomenological parameters of the 
original theory. 

Integrating out the phase field cj) in the action pO| . 
taking account of the fact that i?q=o — -fq=Oi leads to 
the following effective action for vortices 
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(Since the vortex gauge field belongs to the classical sec- 
tor of the theory, the derived effective action contains 
only the = part). The last result, combined with 
the identity p(q)/5(— q) = (q^/Ko)(aq • a_q), yields the 
effective vortex interaction potential: 



where momenta k and q are expressed in units of fcf,. 
For small q (|q| < O.I/ca), numerical evaluation becomes 
troublesome because of the strongly singular character of 
the integrand in Eq. (|14p . However, in the limiting case 
q| — > 0, by replacing dispersion £k with its linearized 
form (ck — > Wh(|k| — fcf,)) we can derive the analytical 
result Pq —> —kf/{6n'^Vb)- Since Vb = |9£k/c^k||^|^j,^, 
the last result can be expressed as 
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Some typical results of numerical evaluation of the re- 
sponse function Pq for |q| > O.lkA are displayed in Fig. [T] 
(where 2kb < k^)- The salient characteristic of these re- 
sults is a knee-like feature at |q| = 2kb, which reflects the 
existence of an effective FS of diameter 2kb. It bears anal- 
ogy to the 2fci?-feature of the paramagnetic spin suscepti- 
bility in 3D, responsible for the RKKY indirect-exchange 
interaction in metals, albeit the 2fcb-feature found here 
comes from the current-current correlator so that it is 
not directly related to the RKKY interaction. The val- 
ues of Pq obtained analytically in |q| limit differ 
just slightly from numerical values at |q| — O.lfcA, indi- 
cating that Pq can be approximated as a constant in this 
numerically-inaccessible region < |q| < O.lfcA. 
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Thus besides the conventional long-range component pro- 
portional to 1/q^, we have an additional component 
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induced by the gapless fermions. 

Response function Pq has to be calculated numerically. 
Yet, prior to numerical computations, Pq can be reduced 
to a two-dimensional principal-value integral 

1 - cos2 e . 



P" 



1,3 



(2^) 
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FIG, 1: Transverse current response function Pq as a function 
of dimensionless momentum, for m = 1,0 and A/|/i| = 2,0, 
Values of kh are given in units of k^- 

By exploiting the rotational invariance of the induced 
potential in momentum space, in real space we obtain 



l^nd(r) = 



1 

2^ 



|q|Vind(|q|)Jo(|qk)rf|q| 



(16) 



where Jo{x) is the zeroth-order Bessel function of the 
first kind. Our numerical calculations for different val- 
ues of relevant parameters (fcb, A) show that the induced 
potential has RKKY-like oscillating behavior, with at- 
tractive character {dV-md/dr > 0) at short distances. To 
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calculate the effective vortex interaction potential, given 
by the sum of T^nd and the repulsive logarithmic part 
Vo(r) = — (Kgns/2mp) ln(fcA'')j we have computed the 
superfluid density by adapting the method of Ref. ^IS^ to 
the special case of a single gapless FS. 

In order to elucidate the realm of validity of our 
effective theory and make contact with experiments, 
it is useful to estimate the physical healing length 
^ = {Snnsam)^^^'^ {o-m being the molecular scattering 
length), with the inverse of the momentum scale k/^. It 
is known that in the strong-coupling BEC regime of a 
superfluid Fermi gas with equal populations of two hy- 
perfine spin components the molecular scattering length 
is given by am — 0.6 a/ (a/ being the scattering length 
between fermionic atoms). For a polarized Fermi 
gas, as shown by Sheehy and Radzihovsky, decreases 
monotonously as a function of the polarization and van- 
ishes at the boundary of first-order phase transition to 
the phase separated state. Therefore, right at the bound- 
ary to phase separation and in the immediate vicinity 
of it the coherence length becomes much greater than 
/c^^, thus making the quantitative implications of our 
theory not directly applicable in this special case. Taking 
o.m{P = 0) in place of a„i, together with typical values of 
Us/n and A/e^ (ef = ^f/(2to), where kp — (Stt^ti)-'^/'^ 
is the momentum scale set by the total fermion density) 
in the BEC regime, we estimate that ^ is of the same 
order as k'^ when |k| ~ 1-100. By taking into account 
the fact that £,/k~^ oc (om/a/)"^/^, we can infer that the 
above estimate is just slightly modified as a result of 
decreasing in the vicinity of the phase separation line. 

As our calculations demonstrate, the effective vortex- 
vortex interaction shows three characteristic types of be- 
havior, i.e. three polarization-dependent regimes. The 
critical polarizations corresponding to the boundaries be- 
tween these different regimes are not universal but de- 
pend on the actual location in the part of the phase dia- 
gram pertaining to the BPl phase. 

In the regime of relatively low polarization, the to- 
tal potential is dominated by the conventional repulsive 
logarithmic part; the effective vortex interaction is repul- 
sive [dVcs/dr < 0) at all distances. The resulting vortex 
phase is accordingly expected to be conventional, with 
triangular vortex arrangement. An example is shown in 
Fig. Hi,. 

In the other extreme - the regime of high polarization, 
the induced potential plays a dominant role at short and 
intermediate distances. This renders the total potential 
attractive at short distances, with pronounced oscillating 
features resembling the RKKY interaction, as illustrated 
in Fig. 

The attractive nature of two-body interaction already 
at short distances suggests an instability of the vortex 
lattice. However, whether this instability really occurs 
is still an open question for the following reasons. The 
physics at distances shorter than the healing length ^ (to 



be discussed in the next section) is not captured by our 
effective theory ; also, the multi-vortex interactions, not 
considered here but certainly allowed as higher orders in 
the effective vortex action, may support unusual vortex 
phases. 

Apart from these two extreme regimes, in a narrow 
window of parameters the total potential is repulsive at 
short distances (r « (2 — 3)fc^^) and becomes attractive 
at intermediate ones (Fig. [3|) . 

Due to the finite range of the RKKY-like induced po- 
tential, the truly long-distance dependence of the effec- 
tive potential (for all polarizations) is governed by the 
infinite-range repulsive logarithmic interaction. 
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FIG. 2: Effective vortex interaction potential in real space 
(in units of A), for m — 1.0 and A/|^| — 2.0. Values of kb in 
both plots are expressed in units of k^: (a) kt — 0.623; 0.724; 
0.846 correspond to polarizations P = 0.155; 0.227; 0.314, 
and (b) kb = 1.500; 1.352; 1.205 correspond to polarizations 
P = 0.763; 0.702; 0.626, respectively. 

In summary, we have found that the vortex interaction 
potential has a remarkable form in the gapless fermion 
superfiuid (BPl), a recently proposed phase for a res- 
onantly interacting imbalanced Fermi gas in the BEC 
regime. The vortex interaction found here contains an 
attractive fermion-induced potential with the RKKY- 
like oscillating character, in addition to the usual re- 
pulsive Coulomb potential. In the conventional super- 
fluids, the point vortices with the repulsive logarithmic 
(2D Coulomb) potential would form the triangular lat- 
tice. [15J Our finding of the novel vortex interaction. 



5 



-0.28 - 




FIG. 3: Effective vortex interaction potential in real space (in 
units of A), for m = 1.0 and A/\fj,\ = 1.0. h = 1.022 (in units 
of /ca) corresponds to P = 0.504. 

therefore, opens up an intriguing question about the 
vortex lattice structure in the spin-imbalanced gapless 
fermion superfluid state. 
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